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ABSTRACT 

We develop a radiation pressure-balanced model for the interstellar medium of high- 
redshift galaxies that describes many facets of galaxy formation at z > 6, including 
star formation rates and distributions and gas accretion onto central black holes. We 
first show that the vertical gravitational force in the disk of such a model is dominated 
by the disk self-gravity supported by the radiation pressure of ionizing starlight on 
gas. Constraining our model to reproduce the UV luminosity function of Lyman-break 
galaxies (LBGs), we limit the available parameter-space to wind mass- loading factors 
1-4 times the canonical value for momentum-driven winds. We then focus our study 
by exploring the effects of different angular momentum transport mechanisms in the 
galactic disk and find that accretion driven by gravitational torques, such as from linear 
spiral waves or non-linear orbit crossings, can build up black hole masses by z = 6 
consistent with the canonical M-a relation with a duty cycle of unity, while accretion 
mediated by a local viscosity such as in an a-disk results in negligible BH accretion. 
Both gravitational torque models produce X-ray emission from active galactic nuclei 
(AGN) in high-redshift LBGs in excess of the estimated contribution fr om high-mass 
X-ray binaries. Using a recent analysis of deep Chandra observations bv lCowie et all 
we can already begin to rule out the most extreme regions of our parameter-space: the 
inflow velocity of gas through the disk must either be less than one percent of the disk 
circular velocity or the X-ray luminosity of the AGN must be substantially obscured. 
Moderately deeper future observations or larger sample sizes will be able to probe 
the more reasonable range of angular momentum transport models and obscuring 
geometries. 

Key words: quasars: general - galaxies: evolution - galaxies: active - galaxies: high- 
redshift - cosmology: theory - X-rays: galaxies 
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1 INTRODUCTION 

Lumi nous quasars have recently been observed out to z > 6 
(e.g., lFanll2006l : iMortlock et all 1201 if ) indicating the pres- 
ence of rapidly accreting supermassive black holes (BHs) 
in the early universe. These active galactic nuclei (AGN) 
are thought to be hosted by extremely rare and massive 
dark matter halos with masses around 10 12 Mq undergoing 
freque nt major merge rs that drive gas toward the galactic 
center |Li et alj|2007r ). 

At the same time, Lyman-break galaxies (LBGs ) have 



been discove r ed a t z > 6 (e.g . . iBouwens et all 120061: 



Bunker et al 
20ld : 



l20ld: 



Bouwens et al 



McLure et all l20ld : iFinkelstein et all 
2011al lbh with typical star formation 



rates (SFRs) of < 1 Mq /yr. An analysis of the resulting UV 
luminosity functions (LFs) shows that these LBGs inhabit 
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~ 1O 1O M halos (|Munodl2012T) that are much more com- 
mon than those that host high-redshift AGN. It is still un- 
certain whether they produce enough high-energy photons 
either to reionize the universe or to kee p it ionized (e.g., 
iMunoz fc LoebhOllljBouwens etld1l2012f ). Yet, for the most 
part, these galaxies are too faint for effici ent high-resolution 
spectra and too small to be fully resolved (|Oesch et al.ll201ol : 
IWvithe fc Loebll201lf ). Thus, observations have, so far, re- 
vealed little about their interstellar media (ISM). 

In the absence of such data, it is often assumed that 
z > 6 LBGs are physically similar to the lower-redshift 
starbursts that have also been found with the Lyman-break 
technique. The fac t that few z ~ 2 LBGs host AGN (e.g., 
ISteidel et al J 120041 ) is often pointed to as evidence that few 
z ~ 6 galaxies sh ould contain supermassive BHs as well (e.g., 
ICowie et al1l2012i 1. Indeed, studi es of X-ray heating in the 
intergalactic medium (IGM) (e.g.. lObl200ll ; lFu"rlanetto et al.l 
2006) have cited the low-redshift AGN LF to argue that 
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AGN are not necessary for cosmic reionization and that 
high-mass X-ray binaries (HMXBs) provide the only sig- 
nificant sources of X-rays. These studies are limited by data 
that may be missing faint AGN in z > 6 galaxies. 

On the other hand, a general theoretical picture 
of galaxy physics is emerging where galactic accretion 
is balanced by star formation, winds, and angular mo- 
ment um transport in marg i nally Toomre-unstable disks 
e.g.. iThompson et al] 120051: iFaucher-Giguere et alJ l201ll: 



Hopkins fe Quataertll201ll ; iHopkins et al.ll2012l ; iDave et all 
and in which gas inflow to ward the center of the galac- 
tic center is a generic predict ion ()Hopkins fc Quataertll2010l . 
l201ll ; iBournaud et all 12011ft . Moreover, recent analysis of 
Chandra observations has begun to place limits on the X- 
ray luminosity of z > 6 LBGs using a stacking method 
|Treister et al.l l201ll : ICowie et af]|2012ft . While the results 
of different studie s have been inconsistent, the conclusions 
by iTreister et al.l ()201ll ) may reflect evidence of ongoing 
and highly obscured black hole growth in faint AGN at the 
centers of these high-redshift LBGs. However, predictions 
for such X-ray observations based on the new paradigms 
of galaxy formation are absent in the literature. With a 
view toward resolving this shortfall, we start from the radia- 
tion pressure-balanced ISM model of lThompson et al . (2005, 
hereafter TQM05) that was shown to describe massive star- 
forming galaxies and AGN at z ~ 2 and extend it to even 
higher redshift by assuming a dust-free ISM and applying 
the appropriate velocity dispersions and metallicities. We 
improve the treatment by considering both the self-gravity 
of the disk and the ejection of material through winds and 
show how the UV mass-to-light ratio and the X-ray lumi- 
nosity are signatures of and place constraints on the galaxy 
fueling rate, wind properties, and angular momentum trans- 
port mechanism. 

In Ej2j we outline the components of our ISM model 
showing explicitly how it differs from previous studies and 
implementations and clearly laying out the free parameters. 
After constraining the wind mass-loading parameter to re- 
produce observations of the UV LF at z > 6 in we explore 
the effects of different angular momentum transport models 
on the central black hole accretion rate in §3] and on the 
resulting AGN luminosity and its obscuration by the dense 
gas in fj5] We present our results in Sj6] where we specifically 
consider the X-ray emission of high-redshift LBGs, compar- 
ing our model AGN outputs with expectations from HMXBs 
and the recent Chandra observations. Finally, we conclude 
with a review and discussion of our conclusions in iJ7| 



2 THE MODEL 
2.1 Overview 

In this subsection, we summarize the basic components of 
the model discussed in the rest of [J2]with fiducial values and 
free parameters. We formulate this model with the goal of 
representing z > 6 LBGs in mind. At lower redshift, it re- 
duces roughly to the model that TQM05 showed accurately 
describes luminous AGN and starbursts at z ~ 2. Here we 
also point out where our formulation differs from that in 
previous studies. 

We start with the basic model of TQM05 with a galac- 
tic disk that connects smoothly to an AGN accretion disk 



around a central BH. We note that, while the disks need not 
be thin, this picture assumes that high-redsift galaxies are 
not spherically amorphous. The size and velocity dispersion 
of the galactic disk are specified by the mass and redshift of 
the host dark matter halo. As described in £|2 . 2 1 gas is ac- 
creted onto the outer edge and, at each radius, may either be 
turned into stars, ejected by winds, or transported toward 
smaller radii. Given their focus on larger galaxies, TQM05 
ignored the loss of gas via winds. However, we include a 
momentum-driven outflow at each radius with a mass loss 
rate proportional to the ratio of SFR to velocity dispersion 
and normalized by the free parameter 770. The timescale for 
the buildup of the disk is long enough that the reservoir 
of existing gas need not be included either as contributing 
to star formation or as a sink for accreting gas. Ultimately, 
the fraction of gas not incorporated into stars or ejected by 
winds is accreted onto the BH. 

The rate of baryonic accretion onto the galaxy ( £12.31) 
and the mechanism for gas transport through the disk ( H2.4I) 
are key elements of the model. We assume that the cold- 
flow model — in which infalling cold gas is not shocked at 
the virial radius but streams freely to the galaxy — gives the 
most reasonable description of the galactic fueling rate in 
z > 6 LBGs and is completely specified by halo mass and 
redshift. Moreover, we consider three possible mechanisms 
for the transport of gas through the disk: (i) an ct-disk model 
with a viscosity parameterized by a constant value a; (ii) 
a gravitational torque-generated infall model induced by a 
linear spiral wave and characterized by an infall velocity that 
reaches a constant fraction, m, of the local sound speed; 
and (iii) a shocked infall model induced by orbit crossings of 
spiral waves in which the infall velocity is a fixed fraction, /?, 
of the orbital velocity. This last process was not considered 
by TQM05 but has be en seen in numerical simulations and 
developed analytically (|Hopkins fc Quataert|[2O10l . l2Ollft . In 
each case, we take fiducial values of a — 0.3, m — 0.2, or 
13 = 0.01. 

Given the global accretion rate, an angular momen- 
tum transport mechanism, and a wind parameter 770, we 
can calculate the density and temperature profiles as well 
as the star formation and infall rates as a function of ra- 
dius in the disk by solving the equation of vertical hy- 
drostatic equilibriurrQ f N2.5[l . in which radiation pressure 
from stars dominates thermal pressure from supernovae and 
maintains marginal Toomre-stability. While we consider the 
effect of additional turbu lent support, eithe r from dense 
clumps in cold streams (iDekel et al.1 2009bl) or disk in- 
stabilities jBurkert et alJ l2010l : iKrumholz fc Burkertl l20ld : 
IBournaud et al J l201ll ). we ultimately ignore both when com- 
pared with pressure from stars ( H2.6|) . Deviating from the 
TQM05 model, we include the self-gravity of the disk in 
the pressure balance and assume a dust-free ISM, taking as 
negligible additional support from pressure on dust grains 



1 Following TQM05, where multiple solutions exist, we select 
the one at the lowest temperature and ignore the multiplicity 
of phases in the ISM that these other solutions might imply. 
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2.2 Balancing Inflows, Outflows, and Star 
Formation 

The TQM05 model assumes a galactic disk that extends to- 
ward small radii and matches onto the accretion disk around 
a central BH. This system is embedded in an isothermal halo 
with velocity dispersion a so that the angular rotation speed 
of the disk as a function of radius is given by 



outflows, we include the effect of a momentum-driven wind 
|Murrav et al.ll2005l ) with 



G M] 



BH 2o2 



(1) 



where Mbh is the mass of the BH. Given the paucity of 
high-redshift data and the uncertainty in extrapolating the 
local M-cr relation to early times, we simply assume the local 
empirical relation 

M BH = 2 x 10 8 a| 00 M (2) 

with (72oo = °7 (200 km/s) (|Tremaine et alj 120021 1 and sub- 
stitute the total halo velocity dispersion for the bulge dis- 
persion. For a 10 10 Mq halo at z = 6, where a w 50 km/s, 
the extra contribution to the rotational velocity from the 
BH is significant only within the central parsec; throughout 
most of the disk, $1 ~ \ /f 2a/r. 

We assume that gas in the galactic disk is optically thick 
to the ionizing UV radiation emitted by stars so that ra- 
diation pressure maintains vertical hydrostatic equilibrium 
against gravity. In the outer regions where pressure from 
stars dominate s, the disk is constrained to be marginally 
Toomre-stable l|Toomrelll964h with 



ttGE e 



= 1 



(3) 



where c s is the total sound speed in the disk, E g = 2 p h 
is the gas surface density, p is the gas density in the plane, 
and the disk scale height is h. Thermal gas pressure and 
radiation from the AGN are allowed to stabilize the inner 
disk to Q > 1. We review this pressure equilibrium in £12.51 
but we refer the reader to TQM05 for more details. 

The model implicitly assumes that the timescale for 
building up the reservoir of gas in the disk is much longer 
than any of the disk crossing timescale, the star formation 
timescale, or the advection timescale. Thus, none of the gas 
accreting onto the disk goes toward increasing the gas mass 
of the disk. For momentum-driven winds, this is true for 
galaxies hosted in halos le ss massive than 10 Mq out to 
z = 10 (jPave et al.ll2012l) . Instead, accreted gas can be 
turned into stars or expelled by winds, as considered by 
iDave et al.li|2012D . but a small fraction may also be accreted 
onto the central BH and power an AGN. This conservation 
is also true on a local level as gas is transported through the 
disk, so that the inflow rate at each radius is: 



M(r) = M dis 



2 7rr'E*(l + 77 wind )dr', (4) 



where Mdisk is the galactic accretion rate of baryons onto the 
outer edge of the disk at a radius of i?disk, E* is the surface 
density of star formation, and 7? w i n d = M w i n d/M* is the wind 
mass-loading parameter. We take -Rdisk to be a fixe d fraction 
X/V2 of the halo virial radius, i? v ir, with A = 0.05 (|Mo et al.l 
1998), co nsistent with the highest-resolution obse rvations 
available (|Oesch et alj l20ld : IWvithe fc Loebll201ll ). While 
TQM05 ignored the reduction in gas accretion rate through 



f?wind = V0 ' 



100 km/s 



IO- 



CS) 



where 770 is a free parameter to which we will return in [J3] 
and e is the efficiency with which stars process matter into 
radiation. In principle, e depends on the IMF and metallicity 
of the stars in the galaxy with an efficiency of eo ~ 10~ 3 for a 
Salpeter IMF from 1-100 Mq at solar m etallicity. However, 
using Starburst99 (|Leitherer et al.l 1 1999T ) . we have checked 
that e differs from eo by only a factor of order unity even 
for a Pop III IMF at very low metallicities. Thus, we will 
assume e = eo throughout the rest of this PaperQ Given this 
efficiency, if the asym ptotic wind velocity at infinity is w 3 a 
(|Murrav et al.ll2005l ). the critical velocity below which winds 
reduce the gas accretion rate more than star formation itself 
is 300 7/0 km/s, implying that most high-redshift galaxies are 
strongly affected by winds. 

On the o ther han d, the semi - analy tic models of 
iRaicevic et al.l ()201 lh and lLacev et al.l (|201ll ) postulate that 
LBGs at z > 6 are the products of mergers, which are prob- 
ably not well-represented by the steady-state model we have 
assumed in equations [T] and [3] At z = 6, the halo merger 
time of 10 10 M R halos is on the order of the Hubble time 
(e.g., iKulkarni fc Loebll2012l ), but the small sizes and high 
densiti es of early galaxies may suppress merger-generated 
effects (Cen 2011). Given the complexities of such a model, 
we leave a detailed analysis of this scenario to future work. 

The galactic accretion rate of baryons onto the disk, 
Afdisk, is a key component of equation [4] since it limits the 
amount of p ossible star forma t ion, w inds, and AGN activity. 
By contrast, iKrumholz et al.l (|2009T ) proposes a star forma- 
tion efficiency set by the properties of molecular clouds. This 
picture can be reconciled with our model if the molecular 
gas fra ction is simply set b y the SFR required to maintain 
Q ~ 1 l|Hopkins et al.ll201ll ). In i)2.3l we consider the galac- 
tic accretion rate due to cold flows. The transport of the gas 
through the disk is also a critical component of the model 
since the transport rate determines whether gas will form 
stars and expel winds (i.e., if the rate is slow) or be ac- 
creted onto a luminous AGN (i.e., if the gas accretes quickly 
through the disk so that the advection timescale is shorter 
than the star formation timescale). H2. 41 considers three dif- 
ferent mechanisms for the transport of angular momentum 
and gas through the disk and their implications. 



2.3 Galactic Gas Accretion 

The accretion rate of baryons onto the disk is an important 
parameter in our model because it sets the amount of gas 
available for star formation and AGN fueling. As we will see 
in i\2A\ the accretion rate also determines the scale height 
required to channel such a mass flow rate through the disk 
at a certain velocity. 

In the cold flow model, baryons stream directly onto 
the galactic disk and are not slowed by an accretion shock 

2 Since only ionizing radiation from stars impacts the dust-free 
ISM in our model, the value of e may be correspondingly reduced. 
However, because of the degeneracy between e and rjo in equation 
\5\ this ultimately has no effect on »? w i n d itself (see footnote [7j- 
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l|Keres et al.ll2005l ; [Pekel et alj|2009a[ h Thus, the accretion 
rate roughly traces the buildup of dark matter in the halo. 
Several studies have investigated fitting formulae for the av- 
erage accretion ra te as a function of only halo mass, Mh a i , 
and redshift (e.g..lNeistein et al.ll200rj ; IMcBride et~aLll2009l ; 



iFaucher-Giguere et al.ll201lh . While all of these studies were 
focused on results at z ~ 2, we assume that the same red- 
shift dependence holds out to z > 6 and gives a reason- 
able approximation of the baryo nic accretion rate at these 
redshifts. Following the work by IMcBride et all l|2009l) . the 
average baryon accretion rate in cold flows at high redshift 
is: 



M disk » 3M /yr 



10 10 M 



1 + z 



h 

0.16 



(6) 

where /b is the cosmic baryon fraction. We note that this 
prescription would be significantly altered in a merger sce- 
nario, where gas may be dumped onto galaxies at much 
highe r rates for shor t periods of time and at chaotic intervals 
fe.g., lLi et alll2007l l. 



2.4 Angular Momentum Transport Through the 
Disk 

The rate at which gas is transported toward the center of 
the disk determines how much is consumed by star formation 
and how much fuels the AGN. We consider three different 
mechanisms that may operate in the disks of high-redshift 
galaxies: a simple Q-disk model, infall mediated by linear 
spiral waves, and non-linear shocked infall due to orbit cross- 
ings. In each case, the mass infall rate as a function of radius 
is given by 



M — 4 7r r hpvin, 



(7) 



where Vi n is the infall velocity of the g function of ra- 

dius. With p determined by equation O the mechanism that 
sets Win also indirectly determines the disk scale height re- 
quired maintain mass conservation while channeling a given 
M through the disk. 

2.4.1 a-disk 

In an Q-disk, the viscosity funneling gas into the center of 
the disk is a local pro cess given by v = a c s h an d described 
by the parameter a l|Shakura fe Sunvaevlll973T ). Following 
TWM05, the infall velocity is 



Win = V 



dlnO 



dr 



(8) 



Unless otherwise noted, we take a fiducial value of a — 0.3. 
However, the inf low can achiev e even lower values of a in 
certain regimes iGammiellioOll ). Since a smaller a would 
only reduce the already miniscule BH accretion rates we 
find for locally viscous disks (see S}, we consider only the 
higher value here. For a thin disk, the infall velocity of an 
a-disk is much less than the sound speed. 



2.4-2 Linear Spiral Wave 

Since a local mechanism for viscosity may have trouble fu- 
eling the central BHs of luminous AGN, TQM05 further 



consider a global-torque model in which the inf all velocity 
reach es a fixed fraction of the sound speed (e.g iGoodmanl 

v in = mc s , (9) 

where m is the Mach number. This type of angular mo- 
mentum transport is produced by gravitational torques 
due to linear spiral waves (LSWs) in pure g as disks 
(e.g.. lKalna!slll97ll: lLv"nden-Bell fc Kalnaisl [l972l see also 
Hopkins fe Qua tacrt 201 1| and references therein). Following 
Goodman! ( 2003) and TQM05, we assume a fiducial value of 
m — 0.2. Thus, the infall velocity in this model is a signifi- 
cant fraction of the sound speed, even for a thin disk. As we 
will see, this allows more gas to accrete onto the BH than 
in the Q-disk case. 



2.4-3 Stellar Torques and Shocked Infall 

Two additional effects will lead to increased angular mo- 
mentum transport and gas infall in galactic disks. First, in 
disks with both gas and stars, wave modes in each com- 
ponent will be of fset by ~ lOdeg due to dissipation in 
the collisional gas dNoguchill 19881 : [Barnes fe Hernquistlll996l : 
iBerentzen et a.l.ll2007f ) . The resulting gravitational torque of 
the stars on the gas drives gas toward the center. However 
this effect is strongly suppressed in systems where the disk 
does not dominate the potential (|Hopk ins fc Quataerd[2O10l . 
l201lf ) . Additionally, non linear wave modes in the g as lead 
to orbit c rossings (e.g., |P apaloizo u fc Pringld Il977l ) , orbit 
trapping (jBinnev fc Tremainel Il987l . e.g..). and ultimately 
shocked dissipation and infall. Both of these processes lead 
to mass infall rates that are independent of the local sound 
speed and, thus, of the lo cal thermodynamics of the disk 
IjHopkins fc Quataertll2011r ). In general, 



(10) 



where j3 is a constant. Because we assumed in H2.2I that 
the disk is a subdominant component of the potential in 
z > 6 galaxies, we do not expect a stellar disk to provide 
a significant torque on the gas. Therefore, we will refer to 
the angular momentum transport mechanism of equation 
1101 as resulting from "shocked infall." On the other hand, 
the formation of disks at the highest redshifts is an open 
problem. Numerical simulations with significant inflow that 
include both effects find t ypical values of /3 ~ 0.001-0.1 
i|Hopk ins fc Quataertll2010h . In what follows, we nominally 
assume (3 — 0.01 and consider results for a range of values 
where appropriate. 

2.5 Pressure Balance 

We assume that hydrostatic equilibrium is maintained ver- 
tically in the disk. For a thin disk, the pressure required to 
balance gravity can be approximated as: 

^2 



2 7T G E g p h + p h? Or 



(11) 



3 It may be the case that BHs are fed by the small frac tion of 
gas with initially low angular momentum l|Goodmanll2003h rather 
than by gas whose angular momentum has been transported out- 
ward. We leave an exploration of this possibility and its incorpo- 
ration into our models for future work. 
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where the first term on the right-hand-side is the contribu- 
tion from the self-gravity of the disk and the second is due 
to the isothermal halo potential. We show that, for constant 
Q, equation 1111 can also be expressed as P gra v = B ph 2 Q 2 , 
where B is a constant. While TQM05 assumes 5 = 1 for 
simplicity, essentially ignoring the disk self-gravity term, we 
show that the contribution of this term is indeed significant. 
Since c 2 — P grav /p, then c s = B 1 ^ 2 hQ. Substituting into 
equation [3] with £ g = 2ph and solving for the density gives 

P=^, (12) 

which is different from the TQM05 derivation by the factor 
of B 1 / 2 . Substituting back into equation 1111 yields P grav = 
(2 3/2 B 1/2 /Q + l)ph 2 Q. 2 = Bph 2 Q 2 from which we find 
(for Q = 1) B = (V2 + V3) 2 w 9.9. While such a large value 
of B implies that the disk self-gravity is an important com- 
ponent of P grav , this is only because the gravitational pull 
of the halo contributes only its small vertical component. 
The potential and total gravitational force throughout most 
of the disk is still dominated by the isothermal halo as in 
equation [T] 

Gravity in the disk is balanced by radiation pressure 
from stars and the thermal pressure of the gas. TQM05 
considered the influence of radiation on grains within a cou- 
pled dust-gas medium. However, in the dust-free ISM that 
we assume for z > 6 LBGs, UV photons impart momen- 
tum directly to the gas through ionizations. In this case, 
the gas surface densities required for optical thickness are 
extremely small (~ 10~ 5 M©/yr). For an ionizing escape 
fraction / csc <C 1, the resulting pressure is of the same or- 
der as that in the single-scattering case of a dusty medium, 
where the dust is optically thick to UV radiation but opti- 
cally thin to the re-radiated IR (e.g.. I Andrews fc Thompson! 
120111 ) . and given by: 

-fstars.rad = £ E* C. (13) 

Note that we have, additionally, ass umed no cancelation of 
oppositely-directed momentum fe.g. JSocrates et alj2008h in 
equation 1 131 

TQM05 showed that the contribution to the pressure 
from a hot ISM phase in supernova bubbles is negligible 
because of the high density of luminous starbursts. How- 
ever, the mechanical deposition of energy from gas swept 
into shells by supernova shocks provides pressure support of 
the same order as Pst ars , ra d. The same argument holds true 
in the comparably high densities of the early galaxies we 
consider here. Thus, we set the total pressure from stars to 
be Pst a rs = 2 P s t a rs, ra d and have checked that this assump- 
tion provides nearly identical results as explicitly including 
mechanical pressure support from supernovae. Because con- 
traction of the disk leads to an increase in star formation 
while expansion leads to a decrease, pressure due to star 
formation is a self-regulating proces^j that maintains Q ~ 1 
in the disk. 

While the disk is dense enough on average to gain mo- 
mentum from every ionizing photon produced, winds or the 

4 Here we implicitly assume the "continuum" approximation dis- 
cussed in footnotes 12 and 14 of TQM05 so that star formation 
feedback continues to maintain Q = 1 even in annuli where the 
SFR is very small and discreteness effects become important. 



AGN may create holes in the gas distribution through which 
radiation can leak and increase / osc . In this case, stabilizing 
the disk will require additional star formation and result in 
less gas available for accretion into the central AGN. How- 
ever, we would expect the resulting leakage to have only a 
small effect on the mechanical pressure supplied by super- 
novae. Therefore, holes should not change P s t a rs by more 
than a factor of two. The same argument can be applied to 
using more precise values of e for pressure generated exclu- 
sively by ionizing radiation. 
Thermal gas pressure, 



pk B T 



(14) 



also provides support against gravity according to its density 
and temperature. Here, m p is the proton mass, and, in the 
star-forming region of the disk 



(15) 



so that the mid-plane temperature, T, is approximately the 
effective temperature of the radiation^ We also ignore the 
effect of CMB heating on the gas. The CMB temperature at 
z = 6 is about 19 K, which can be higher than our modeled 
gas temperature at some radii and in some regions of param- 
eter space. However, the Compton heating time for the CMB 
is long, and formally, there are no metal lines through which 
it can interact with the gas in our dust-free model. Still, the 
CMB may have an important effect on the temperature of 
high-redshift molecular clouds and certa inly on our ability 
to ob serve radio lines from these objects |Obreschkow et all 
l200Sh . 

Vertical hydrostatic equilibrium implies 

P g r a v — Pst a rs + P gas . (16) 

As in the TQM05 treatment, once P gas comes to dominate, 
star formation is no longer required to support the disk and 
becomes negligible. As gas accretes to even lower radii, it 
loses the gravitational energy that powers the AGN. How- 
ever, since star formation can no longer deplete the supply 
of gas to the central BH, we ignore the question of how and 
whether this extra energy heats or supports the disk. 



2.6 Turbulence from Dense Gas Clumps and Disk 
Instabilities 

In this subsection, we consider turbulent pressure sup- 
port due to t he impact of dense clumps in a cold flow 
on the d isk (IDekel et alj l2009b|) or instabilities within 
the disk dBurkert et all l20ld ; iKrumholz fc Burkertl l20ld ; 
iBournaud et al.ll201ll ). which have been postulated to be 
the dominant mechanisms maintaining Q = 1. 

Since the angular momentum of cold flows i s thought to 
be responsible for the bui ldup of rotating disks (IDekel et al.l 
l2009ri [Stewart et al.|[201ll ). we assume that dense clumps 



5 We note that our d ust-free formulation differ s from the treat- 
ment in TQM05 and ISirko fc Goodmanl ||2003|) where the mid- 
plane temperature T 4 = T 4 ff (3r/4 + r~ 1 /2 + 1), where r is the 
IR optical depth. Additionally, a lower value of e would not sig- 
nificantly affect the gas temperature because the corresponding 
reduction in pressure would be raised to the 1/4 power. 
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within the streams impact and deposit their kinetic energy 
as turbulence into the disk at its edge. This turbulence is 
transported inward over the advection timescale ~ Rdisk/fm 
but decays over the disk crossing time ~ h/c B . For an LSW 
model of gas transport (Eq. [9J, the decay time is a frac- 
tion ~ mh/ Rdisk of the advection time and even smaller 
in an a-disk model. In a shocked infall model, this ratio is 
~ /3ah/(cs Tidisk)- Thus, for a thin disk, m <C 1, or f3 <C 1, 
cold clump turbulence decays quickly and does not affect 
the disk beyond its outermost rim. 

A related source of turbulent disk support may come 
from the energy released by the gas as it spirals inward 
through the disk and down the potential gradient. If all of 
this energy is quickly radiated, then pressure of this magni- 
tude is sub-dominant compared to stars, except very close to 
the galactic center. Otherwise, if the turbulence remains un- 
thermalized, it may dominate and stabilize the disk as sev- 
eral studies that ignore pressure from stars claim in the lit- 
erature fe.g.lBurkert et al ] |2010t iKrumholz fe BurkertfeOld ; 



iBournaud et al. 2011). Formally, this pressure may suppress 
all star formation in our model, but such an occurrence de- 
pends sensitively on, e.g., the details of the turbulent dissipa- 
tion. However, dominant instability-driven turbulence is in 
tension with other works that show star formation as partic- 
ipating directly in disk support and being re gulated by the 
amount of cold inflow into the galaxy (e.g. iHopkins et al.l 
l201ll : lDave et alJl2012h . Our goal is not to definitively settle 
this question, which must be left for numerical simulations, 
but since our work implicitly assumes the latter theoretical 
framework, we here ignore the contribution of lost gravita- 
tional energy to the turbulent support of the disk. 

Finally, we note that additional turbulence in a major 
merger scenario-not considered here-may maintain Q > 1 
and allow gas to bypass star for mation and fu nnel rapidly 
into a central, bright AGN (e.g.. lLi et 1111120071) . 

2.7 Dust Content and Metallicity 

Our treatment has assumed that the ISM of z > 6 galax- 
ies are dust-fre«B However, the dust content and metallicity 
of high redshift galaxies is not yet certain and widely de- 
bated. While metals and dust have clearly be en observed 
in z ~ 6 quasars (e.g. JWang et al.ll2007l . l2011bl lah. the blue 
spect ral slopes of LBGs i ndicate t hat they may bedust- 
free l)Bouwens et al.l l2010l . but see iDunlop et all l|2012h ). 
Theoretically, s ome studies postu l ate metal-free population 
HI stars (e.g . . iTrenti fc; Stiayellil 120091: ITrenti et all 120091 : 
|johnsonl201ol ; ICai et al.ll201ll;IZackrisson et alj|201lh or top- 



heavy IMFs (e. g., iRaicevic et all l201ll ; lLacev et all l201ll ; 
iKim et aT1l201ll ) in galaxies at z > 6. Others find that nor- 
mal stellar populations are sufficient to reproduce the mea- 



There may, if fact, be significant amounts of dust confined to the 
cores of giant molecular clouds where star formation takes place 
and on which radiation press ure can act to supp ort the clouds 
against gravitational collapse llMurrav et al ]|201of) . However, this 
does not contradict our dust- free assumption, which applies to the 
large-scale regime where we have taken the gas and star formation 
to be smooth and homogeneous at a given radius. A lack of dust- 
mixing may be a natural consequence of the density of high-redsift 
galaxies into which supernova bubbles have difficulty expanding 
(see g22}. 



sured UV LF but disagree about the amount of dust extin c- 
tion required (Salvatcrra et al.ll201ll ; iFinlator et alll201ll ). 

For our purposes, whether momentum is supplied to the 
gas directly from ionizing radiation or via UV incident on 
dust grains, the resulting radiation pressure will not be sig- 
nificantly different from that given in equation [T3] as long as 
the dust is optically thin to the re-radiated IR. In luminous 
AGN and starbursts, this reprocessed radiation can pro- 
vide addition al pressure support for m aintaining hydrostatic 
equilibrium (|Sirko fc Goodmanl liool , TQM05), increasing 
Rstars.rad by a factor of (r + 1), where r = cE E /2 is the 
optical depth of dust to IR radiation and n is the Rosseland 
mean opacity in the optically thick limit. Since k gives the 
absorption cross-section per unit gas mass, it is ultimately 
k that reflects the assumed metallicity and dust-to-gas ra- 
tio of the system. Despite the debate surrounding the exact 
dust content of z > 6 galaxies, it is likely that they have 
significantly lower dust opacities than the ones appropriate 
for protoplanetary d isks in the Milky Way (|Bell fc Linlll994l ; 
ISemenov et~aLll2005t l and used in TQM05. 

Our dust-free approximation is further justified by the 
small velocity dispersions and infall rates of interest. For an- 
gular momentum transport mediated by LSWs, combining 
equations [7] and [9] yields 



S e = 



oM 



7t 2 GQmr 2 



(17) 



Assuming the simple, low-temperature opacity model used 
by TQM05 where k — kq T<j, with dust temperature Tij and 
1 cm 2 g^ 



ko « 2.4 x 10 4 cm 2 g 1 K 2 , the resulting optical depth is 



4.2 



300 km/s 



M 



320 M Q /yr 



200 pc 



100 K 

(18) 

While the luminous starbursts and AGN considered by 
TQM05 have o ~ 300 km/s and M ~ 3OOM /yr at the 
edge of the disk, z = 6 LBGs are significantly smaller sys- 
tems with a ~ 50 km/s and a cold-flow accretion rate of 
M ~ 3 Mq /yr. This results in a much lower IR optical 
depth, even for similar opacities. Repeating the argument 
for shocked infall and scaling to values appropriate for high- 
redshift LBGs, the optical depth is 

"" S ' 50 km/s 



M 

3M /yr ) V200pc 



f T d 



VlOOK 
(19) 

While, in this case, r is of order unity at the edge of the 
disk, we will see in [|4]that M/r drops rapidly just inside 
the edge, so that t< 1. 

The uncertainty in the dust content of z > 6 galaxies 
is also an issue for determining the relationship between the 
intrinsic UV radiation (at 1500 A) from star formation and 
that observed. In we consider two possibilities for the 
amount of UV dust extinction present in high-z galaxies as 
we constrain the wind parameter 770. However, we antici- 
pate that upcoming observations with the Atacama Large 
Millimeter Array (ALMA) will shed more light on this point. 



3 CONSTRAINTS FROM THE UV LF 

The model in [J2]sets the average star formation and AGN ac- 
cretion rates for galaxies hosted in halos of a given mass and 
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redshift. In this section, we use current measurements of the 
UV LF to (1) determine the appropriate halo masses corre- 
sponding to 2 = 6 LBGs and (2) constrain the wind param- 
eter 770. This will ensure that our model galaxies will have 
the abundances and total star formation rates observed. 



3.1 Average M/L 

iMunoz fc Loebl (|201lh and lMunod (|2012l ) inferred the evolu- 
tion in the mass-to-light ratio, M/L, based on UV LF obser- 
vations. In their model, each galaxy has a distribution of pos- 
sible luminosities with probabilities depending only on the 
mass of the host halo given by dP(Li5oo|Mh a io)/dlogI/i5oo- 
Here, we invert this relation using the Bayesian method, to 
obtain the probability distribution of halo masses associated 
with a specific luminosity at 1500 A, L1500: 

dP(M ha i |ii 5 oo) , dn(M halo ) dP(Li5oo|M ha io) 

OCeAF(Mhalo) 



dlogMhalo 



dlogI/1500 



(20) 

Here, dn/dMhaio is the halo mass function, and eaf is the 
fraction of halos that host galaxies and depends on a criti- 
cal suppression mass, Af sup p. Since larger halos are built up 
from smaller ones, eaf < 1 even for some halo masses larger 
than M supp . However, if the minimum halo mass capable of 
hosting a galaxy is Af supp = 10 8 Mq , then eaf ~ 1 for all 
currently observed z ~ 6 LB Gs. Fitting the mode l to the 
UV LF of i-dropouts at z ~ 6 (|Bouwens et al.|[2007l ) reveals 
that dP(I/i5oo|A^haio)/dlogLi50o is a log-normal distribution 
with mean given by 



(logLisoo) w 27.19 



10 10 Mr. 



(21) 



at z = 6 and roughly constant stand ard deviation ol ~ 0.25 
jMunoz fc Loebll201ll ; lMunoj|2012l ). 

Using this model, the average halo mass corresponding 
to the faintest HUD F i-dropout with an a bsolute magnitude 
of Mab = -17.3 dBouwens et all 120071 ) is approximately 
10 10 Mq with a distribution in log-space having a standard 
deviation of a m ~ 0.24. At z = 6, this average mass corre- 
sponds to a velocity dispersion of a ~ 46km/s. The mass 
of the halo hosting each galaxy specifies the average baryon 
accretion rate for a given accretion model. It also sets the 
black hole mass based on our assumed M-a relation (Eq.[2]); 
a 10 10 M Q halo has a central black hole of » 5.6 x 10 5 M . 
While these intermediate mass BHs are still theoretical, they 
can be built up by z = 6 from 100 Mq seeds at z = 20 with 
a constant Eddington ratio of approximately 0.5. 



the discussion in lMufioz fc Loebl l201ll ). however dust may 
diminish this intrinsic UV brightness. To match the average 
observed luminosity of a 10 Mq halo required by equation 
[21] (~ 1.6 x 10 27 ergs/s/Hz), we consider two possibilities: 
(1) that the amount of dust extinc tion is w 0.18 dex as 
determined bv lBouwens et ail l|2007t ) from measurements of 
the UV continuum slope and (2) that rjo = 1, where the 
difference between the resulting UV LF and that observed 
is a consequence of a flexible amount of dust extinction. 
In both cases, we ignore the slight inconsistency between 
the presence of dust extinction and our dust-free pressure 

balance model from ^2] 

The first case was considered by Munoz (2012). Since 
the amount of dust extinction at z as 6 is fixed, then for a 
given Mdisk, Vo is constrained to provide the required aver- 
age luminosity in a 10 10 M Q halo. For the cold-flow accre- 
tion rate in equation [6] the result is rjo ~ 4, where we have 
ignored the remaining slight mass dependence. This is con- 
sistent with results from numerical simulations that require 
superwinds with rjo — 3 (in our notation) to mat ch observa- 
tions at all redshifts in the WMAP 3 cosmology l|Dave et al.l 
2006; lOppenh eimcr fc Pavel l2008h . If the average mass ac- 
cretion rate were significantly larger, as may be the case 
in a scenario where z w 6 LBGs are identified with ma- 
jor mergers, the wind mass-loading factor would need to be 
correspondingly higher to produce the same average UV lu- 
minosity. 

In the 'flexible dust' case, we set r]o = 1 and cal- 
culate the dust extinction coefficient required to produce 
the correct observed luminosity. Given this wind parameter 
at z w 6 in a cold-flow accretion model, a 10 10 Mq halo 
must host galaxies that generate an average SFR of about 
1 Mg/yr. This implies ~ 0.68 dex of dust absorption to pro- 
duce the correct avera ge observe d UV luminosity inferred 
from the UV LF bv lMunozl l|2012l ). Again, we note that the 
additional accretion rate in a merger scenario could be com- 
pensated for by a significantly increased amount of dust. 

The true physics involved in the formation and obser- 
vation of high-redshift galaxies likely involves a combination 
of strong winds and some dust extinction. While we adopt 
a fiducial value of rjo = 4 in the rest of this Paper, we will 
also frequently show results for rjo = 1 given the uncer- 
tainties in calculating the dust extinction from UV contin- 
uum slopes. These values of 770 correspond to critical wind 
velocities-below which winds remove more gas from the disk 
than does star formation-of 1200 km/s and 300km/s, re- 
spectively. This implies mass-loading factors for 10 10 M 
halos at z — 6 of r/wind ~ 8 and 2, respectively^] 



3.2 Constraining Winds 

We want to consider parameter sets of our model that rea- 
sonably reprod uce the M/L a nd its evolution with redshift 
determined by iMunozl (|2012h from fits to the UV LF at 
z > 6. Since we expect the BH growth rate to be small 
compared to the galactic accretion rate, equation[4]becomes 

Mdisk « (1 + ?7wi„d) Mi, (22) 

so that the SFR is set by the galactic accretion rate and the 
wind mass- loading parameter is given by equation [S] At the 
same time, a Salpeter IMF gene rates a luminosity per SFR 
of 8 x 10 27 ergs/s/Hz/(M Q /yr) (|Madau et al.lll99Sl , but see 



4 THE BLACK HOLE ACCRETION RATE 

As gas funnels toward the center of the disk, most is con- 
sumed by star formation or ejected by winds. What remains 

7 Although we have expressed the procedure in this section as 
fitting the parameter rjo for fixed e, equation [22] shows that ?? w i n d 
is really the more fundamental parameter. Thus, as mentioned in 
footnote [2] using a more realistic value of e may result in different 
values rjo because of the degeneracy between the two in equation 
[5] However, this will ultimately have no effect on the gas dynamics 
in our model. 
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Figure 1. The gas accretion rate (thin) and SFR (nr 2 S*; thick) 
as functions of radius in a galaxy hosted by a 10 10 Mq halo at 
2 = 6 with rjo = 4. Short-dashed (red), solid (black), and long- 
dashed (blue) lines correspond to correspond to a-disk, LSW, 
and shocked infall models with a = 0.3, m = 0.2, and = 0.01, 
respectively. 



as the fraction of the age of the universe the AGN must be 
accreting at a constant rate to produce a BH of the given 
size. The highest growth rate of the models we consider is 
produced by shocked infall with /3 = 0.1 and 770 = 1 and re- 
quires a low duty cycle of approximately 0.001. The effect of 
variations in m and /3 for the LSW and shocked infall mod- 
els, respectively, as well as changes in wind strength on the 
rate at which gas is transported through different parts of 
the disk can be seen in Figure [21 On the other hand, a-disk 
models cannot come close to building up such a black hole 
without assistance from an additional phase of intense feed- 
ing. Given the insignificant BH accretion rate produced by 
the local viscosity of an a-disk, we consider only LSW and 
shocked infall models in subsequent sections of this Paper. 

The effects of varying the host halo mass or redshift on 
Mbh are complex since they depend on the details of the 
transition to the AGN accretion disk. However, numerical 
results in the LSW and shocked infall cases are shown in 
the top panels of Figure [3] While, the cold-flow galactic ac- 
cretion rate at a given redshift increases roughly as Mh a i , 
the AGN accretion rate is sublinear with respect to Mhaio 
in LSW and shocked infall models. 

Finally, we note that our calculations ignore the effect 
of strong gravitational wave recoil during BH mergers. In a 
scenario where high-redshift LBGs are the result of mergers, 
the newly merged BH may be kicked away from the dense 
center of the accret ion disk and its growth rate correspond- 
ingly reduced (e.g.. Blecha fc Loeblfiooil ; iTanaka fc Haimanl 
120091: iBlecha et al.ll201ll ). 



is accreted onto the nuclear BH. The rate at which gas is 
transported through the disk critically determines how much 
is available to feed the AGN. The more rapid the trans- 
port, the more gas bypasses star formation and expulsion by 
winds. We begin by numerically solving equationUand then 
investigate the qualitative behavior of approximate analytic 
solutions for each angular momentum transport model. 

Our numerical results show that, as expected, stronger 
winds and slower mass accretion velocities lead to steeper 
declines in the amount of gas transported to smaller radii. 
If the decline is too steep (as in the cv-disk case), star for- 
mation uses up nearly all of the gas at the very outer edge 
of the disk resulting in a ring-like morphology in the UV. 
Conversely, a more even star formation profile is achieved if 
the gas transport is rapid. These results are shown in Figure 
[1] As the density of the disk increases towards its center, gas 
pressure begins to dominate stellar pressure. No new stars 
form inside this radius, and the supply of gas inflow remains 
constant for r < i? agn so that Mbh ~ Af(i? agn ). In Figure [T] 
this occurs where the SFR turns down steeply at very small 
radii. For m = 0.2 in an LSW model with 770 = 4, J? ag n is 
on the order of a few parsecs. 

As depicted in Figure [TJ a-disk models, in which an- 
gular momentum transport is mediated by local viscosity, 
supply more than five orders-of-magnitude less gas to the 
AGN than LSW or shocked infall models for typical val- 
ues of a, m, and /3. We speculate that this dichotomy may 
represent two distinct modes of black hole growth. Figure 
[2] shows that, for Af ha i = 10 10 M o , an LSW model with 
m = 0.2 and 770 = 1 can build up the 5.6 x 10 5 M Q BH 
given by our assumed M-a relation (Eq. [2]) by z ~ 6 with 
a duty-cycle of about unity, where we define the duty cycle 



4.1 Linear Spiral Wave 

Exploring the behavior of the solution to equation [3] under 
some simplifying assumptions can give us qualitative insight 
into the influence of different angular momentum transport 
mechanisms on the BH accretion rate. In the regime where 
fst a rs 3> fpg a s, the vertical hydrostatic equilibrium (Eq. ll6|) 
yields 



Bph 2 si 2 
27c 



(23) 



with p given by equation 1 121 Additionally, combining equa- 
tion [7] with an expression for the mass infall velocity gives 
the disk scale height as a function of radius that is required 
to funnel gas inward at a rate M. For an LSW model with an 
infall velocity given by equation^ we find that h/r oc M 1 ' 2 . 
Further substituting back into equation [4] yields 



M(r) w M d 



A 2 



— dr', 
r 



(24) 



where 



A 2 



B 1 ' 2 o(l + 



?7wind J 



8.3 



2 v / 2ecm 
a 

50 km/s 



l + ?7» 



(HL)- 1 
\0.2J 



(25) 



Equation [24] is an integral equation for M whose solution is 



M(r) w M disk 



^disk 



(26) 
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100 200 300 100 200 300 

r/pc r/pc 

Figure 2. The gas accretion rate as a function of radius in a galaxy hosted by a 10 10 Mq halo at z = 6. Thin lines assume tjq = 1, 
while thick lines assume rjo = 4. Solid curves show numerical results, while dashed curves show analytic solutions from Eg. 1261 and 1321 in 
the regime where pressure from stars dominates the hydrostatic equilibrium. The left panel shows results for LSW infall with m = 0.2 
(black) and m = 1 (red), while the right panel assumes shocked infall with blue, green, and cyan curves denoting = 0.01, 0.1, and 
0.001, respectively. 



This solution was derived in Appendix D of TQM05, where 
their value of A2 differs from ours because of our inclusion 
of winds and the disk self-gravity. 

Equation [26] reveals how the input parameters of our 
model influence the general behavior of the disk. The infall 
rate is such that gas is available for star formation through- 
out the disk, as shown in Figure [1] The M drops off less 
steeply and the BH grows more rapidly for smaller values of 
A2. As we found numerically, this is achieved for smaller val- 
ues of r/o (i.e., less expulsion by winds) and larger values of m 
(i.e., faster infall). Moreover, including the self-gravity of the 
disk also results in slower BH growth since more gas must 
be turned into stars to maintain hydrostatic equilibrium, 
while the stronger supernovae increase the BH growth rate 
by allowing a lower SFR to produce the same hydrostatic 
pressure. The dependences of the BH growth rate on halo 
mass and redshift are complicated since they enter through 
Mdisk, -Rdisk, and i? agn , as well as the precise value of A2. 



We perform the above procedure in i]4.2l and £|4.3l for a- 
disk and shocked infall models of angular momentum trans- 
port. 



4.2 a-disk 

For an a-disk model in the regime where pressure from stars 
dominates the hydrostatic equilibrium, equations [7] and [8] 
combine to give h/r oc M 1 ^ 2 , and equation [4] becomes 



M{r) » M dis 



At 



M 2/s 



dr' 



(27) 



where 
Ai = 



B 5 ' a a 2 (1 + T7 wind ) 
V2ecG 1 / 3 Q 1 / 3 a 

1/3 / o 



74M 1/3 



yr 



50 km/s 



and r/ w i n d is given by equation [5] Differentiating both sides 
of equation [27] with respect to r and solving the resulting 
separable differential equation for M, we find 



M(r) « M dis 



Ai In (Rmsk/r) 



3M, 



>V3 



(29) 



Because of the cubic power on the right-hand-side of equa- 
tion [29] the value of M drops quickly when r < J?disk, as 
seen in Figure [1] The result is a ring-like galaxy morphol- 
ogy where all of the star formation takes place in the very 
outer parts of the disk in contrast to the smoother star for- 
mation distribution of the LSW case. As r — > in equation 
1291 M — > —00. Of course, at some radius, the pressure from 
stars no longer dominates, and the assumptions under which 
we derived equation [57] break down. Indeed, M ^ at ev- 
ery radius in our model. However, the unbounded nature of 
equation [29] highlights the rapid decline of M with decreas- 
ing radius. Larger values of Ai result in an even steeper 
decline. Consequently, the expulsion of more gas via winds 
(i.e., larger rjo) results in slower BH growth, while a faster 
momentum transport rate (i.e., higher a) allows gas to by- 
pass star formation and winds and increases the BH accre- 
tion rate, as expected. 
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4.3 Shocked Infall 

Finally, in a shocked infall model, h/r ex. M, and 

^disk jf^~ 2 



Af (r) « M dis 



■ dr', 



(30) 



where 



A 3 



B 1/2 GQ(1 + 



^?wind J 



16^2eca 2 /? 2 
7O.M0 _1 yr 



v 50 km/s 
and the solution to equation 1301 is 

M disk 



1 + 7?win 

9 



0.01 



'(31) 



1 + A s Mdisk In (7?disk/r) 



(32) 



Because of the large range of potential values of j3 in shocked 
infall, the availability of gas to form stars inside the outer 
edge of the disk depends on the value of A3. When A3 ^> 
100, the amount of available gas falls off steeply with de- 
creasing radius (though, it's not as sudden a drop as in the 
a-disk case) and a ring of star formation develops at the 
outer edge. On the other hand, the more rapid infall rate 
of a disk with A3 <C 100 leads to a smooth, potentially flat 
SFR profile and a higher BH growth rate. As we have already 
seen, BH growth is accelerated by an increase in the infall 
rate (i.e., higher values of 0) or a decrease in the expulsion 
rate from winds (i.e., smaller r]o). 



5 THE AGN LUMINOSITY 

We now turn our attention to the luminosity generated by 
the BH accretion discussed in We consider the BH ra- 
diative efficiency and bolometric X-ray correction (3571} as 
well as obscuration along the line of sight ( 35.20 . 



o 



60 
O 




O 



8 9 10 11 128 9 10 11 12 
log M hal0 /M log M halo /M 

Figure 3. The AGN accretion rate (top panels) and Eddington 
ratio (bottom panels) generated by the BHs in our model at 2 = 6 
as functions of host halo mass. Thick and thin lines assume T70 = 4 
and 1, respectively. Solid (black) and dotted (red) lines in the 
left-hand panels show results for LSW models with m = 0.2 and 
1, respectively, while long-dashed (blue), dot-dashed (green), and 
short-dashed (cyan) curves denote shocked infall models with /3 = 
0.01, 0.1, and 0.001, respectively. For reference, the thick magenta 
line marks Mhalo °c M agn , while the thin magenta line shows the 
accretion rate required to build up a black hole commensurate 
with our assumed M-a relation by z = 6 with a duty cycle of 
unity. 



5.1 Intrinsic X-rays 

The lost gravitational energy of the infalling gas in the inner 
accretion disk powers the central AGN so that the bolomet- 
ric flux emitted at each radius is 

Wr) = ^fl-^) (33) 
Sty V V r I 

If we take the radius of the inner edge of the disk to be 
that of the innermost stable circular orbit (ISCO) around 
a Schwartzchild black hole, Ri n — 6GMbh/c 2 , then the 
radiative efficiency of the AGN, e agn = L agn /(M agn c 2 ), is 
a constant e agn « 8.33% independent of any other model 
parameter[f| For a rotating central black hole, the radiative 
efficiency may be as high as tens of percent. 

Using this efficiency factor, we calculate the bolomet- 
ric AGN luminosity, L agn and the Eddington ratios shown 
in the bottom panels of Figure While models with very 
fast angular momentum transport seem to produce super- 
Eddington radiation, especially in small halos where the cen- 
tral BHs are smaller, this situation is probably unphysical. 

8 This result ignores a general relativistic correction of order 
unity to the radiative efficiency. 



These ratios could be achieved if the emission were colli- 
mated into jets, but it is more likely that the radiative ef- 
ficiency would drop below that calculated above as outgo- 
ing photons are trapped by the infalling gas and ultimately 
dragged into the BH or that the radiation pressure would 
drive an AGN wind. 

While the bolometric AGN luminosity is proportional to 
the AGN accretion rate given a constant radiative efficiency, 
X-ray luminosity is generated by Compton scattering pro- 
cesses that are difficult to predict from first principles. The 
calculation becomes more tractable if we assume a typical 
AGN spectrum to determine the X-ray flux from the bolo- 
metric luminosity. While a handful of quasars have bee n 
discovered at z > 6 (e.g.. lFanll2006l : LMortlock et al.l 1201 lh . 
these objects are quite rare. Since substantial populations 
have not been built up to develop reliable template spec- 
tra at high redshift, we will instead deviate from our prin- 
ciple of not calibrating based on local observations and as- 
sume the bolometric correction for X-ray luminosity given by 
iMarconi et"aH (|2004h . Where necessary, we adjust the given 
correction values of a particular band-for either a slightly 
different energy range or a band corresponding to the rest- 
frame at a higher redshift -by assuming an X-ray spectral 
slope of 1.9. 
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5.2 Obscuration 

Because high-redshift galaxies are much denser than their 
low-redshift counterparts, the X-ray luminosities may be 
absorbed by high column densities of hydrogen. Here we 
consider the amount of obscuration in terms of the fraction 
of impeded sight lines and the amount of gas available. To 
give an approximate description of the central accretion disk 
properties, we assume that vertical hydrostatic equilibrium 
is maintained by radiation pressure in the single-scattering 
limit (see TWM05) generated by gravitational energy lost 
via accretion and no longer constrain Q = 1. 

If X-rays originate from a radius of order the inner edge 
of the disk, _Ri n , then the fraction of obscured sight lines is 
approximately given by h- ln /Ri n , where h ln is the disk scale 
height at i£j n . While a variety of model parameter sets can 
lead to geometric ally thick disks in the TQM05 formalism 
l|Ballantvnell2008l 1. we find that all of our model disks are ex- 
tremely thin in their centers. However, even in a model with 
hi n /Ri n <C 1, sufficiently high X-ray luminosities from the 
centr al source can gene rate a obscuring wind from a dusty 
disk (|Chang et al.ll2007l ). Despite the fact that we have em- 
ployed a dust-free framework for our calculations, our SFRs 
and BH accretion rates would be unchanged if we added a 
significant amount of dust to the inner BH accretion disk. 
Given the contrast bet ween the metals and dust o bserved in 
all z ~ 6 quasars (e.g.. I Wang et al.ll2007l . l2011bllah and the 
blue spectral slopes of the stellar continua (jBouwens et all 
2010), this two-zone dust model may indeed be appropri- 
ate, though of course, the galaxies hosting observed quasars 
are much more massive than the LBGs considered here. In- 
deed, even with a dusty central disk, the intrinsic X-ray lu- 
minosities generated in our models are well below the level 
required to produc e an inflated, unbound disk photosphere 
ijChang et alfeOOTl ). Thus, we do not expect this mechanism 
to significantly increase the fraction of obscured sight lines. 

Nevertheless, to give a sense of its potential effect, we 
consider a generous, toy-model estimate of X-ray obscura- 
tion by neutral hydrogen. In this model, we assume that 
X-ray production proportional to the locally emitted grav- 
itational energy is distributed throughout the plane of the 
BH accretion disk and must escape up through a local disk 
scale height, h. At each radius, the X-rays are spectrally 
filtered by an energy-dependent factor of e ~ s<TH /( 2m p> 5 
where oh = 10~ 17 cm 2 (E/13.6 eV) -3 is the hydrogen photo- 
ionization cros s -secti on. While the bolometric correction of 
iMarconi et"aH (|2004h is calculated based on the total bolo- 
metric emission and, therefore, does not apply to the emis- 
sion at a particular radius, we set the radially-dependent 
bolometric correction to be the same in both the extincted 
and un-extincted cases and assume a fixed cross-section for 
each rest-frame band calculated from the central band en- 
ergy, E. The resulting effective column density is 



N(E) 



In 



Jo 



Eg rav v dv 



a K (E) e- E s -H(S)/(2m p) Fgr^rdr 



(34) 

where the intrinsic emission in a band with average energy E 
is suppressed by a factor of e~ Nf - E ^ °" H ( B ) . We find that a sig- 
nificant amount of X-ray emission generated in the plane of 
the AGN disk is extincted by the dense gas. Only a fraction 
~ 10~ 4 -0.1 of the emitted X-rays in the 3.5-14 keV band 



escape to the observer. The precise amount depends on the 
model parameters, and the exponential dependence on col- 
umn density makes the range large. In general, higher values 
of m or /), while generating larger BH accretion rates and 
Eddington ratios (see Fig. [3]), also result in larger column 
densities. We note that adopting 770 = 1 affects the column 
density only minimally in comparison to changes in the an- 
gular momentum transport rate. We note that the effective 
column density given in equation [34] is significantly higher 
than the average surface density of gas within the accretion 
disk, i.e., E/(2m p ). The difference between E/(2m p ) and 
N(E) is analogous to the difference between the average and 
effective optical depths in studies of the Lya forest. If the 
AGN accretion disk were uniform in density and brightness, 
these two quantities would be equal. As it is, the increase 
in -Fgrav and S toward the disk center also make the results 
independent of the precise value of R asn - 

In the following section, we compare this highly ob- 
scured case to the case of unobscured emission appropriate 
if production of X-rays is near R[ n and the disk is either 
geometrically thin or viewed face-on. 



6 RESULTS 

In the previous sections, we built up and calibrated a model 
that allows us to calculate galaxy disk properties, SFRs, and 
X-ray luminosities as functions of host halo mass and red- 
shift for choices of the angular momentum transport mech- 
anism and wind parameters. In this section, we compare 
the resulting X-ray emission from z > 6 LBGs with that 
produced by star formation ( £16, 1[) and with Chandra obser- 
vations (3631. 



6.1 Comparison with HMXBs 

We first compare the X-ray luminosity from AGN in z — 6 
LBGs to that owing to HMXBs, the dominant mechanism 
for X-ray productio n in star-forming ga laxies without AGN 
at low redshift (e.g.. iGrimm et alJl2003 ). Since HMXBs are 
short-lived, this mechanism is directly tied to the galactic 
SFR able to replenish them, a quantity straightforwardly 
calculated in our model. The effect is usually parameterized 
by the ratio of X-ray luminosity to SFR: cx = Lx/SFR. 
Measurements of the X-ray background and of individual 
galaxies at low redshift find a reasonable range for cx such 
that the distribution of possible values is log-normal with 
a m ean of (logcx) = 39.4 and a standa rd deviation of 
0.4 (|Mineo et alj|2012l ; iDiikstra et afl|2012h . If the AGN in 
LBGs at z > 6 produce a significant amount of X-rays per 
SFR compared to HMXBs, it could affect the heating and 
ionization-state of the I GM and the 21 cm radiation sign al 
during reionization (e.g.. [Oh|[200ll ; iFurlanetto et al.ll2006l ). 

Figure [4] shows the ratio of X-ray luminosity to SFR 
generated by the AGN in our models compared to the es- 
timated range of cx for HMXB consistent with low-redshift 
observations. We have taken Lx in the 0.5-8 keV (rest- 
frame) band, which we obtained by combining results us- 
ing bolometric corrections for the 0.5-2 keV and 2-8 keV 
bands. We find that most models produce intrinsic X-rays 
from AGN in excess of that from HMXBs except in the 
/3 — 0.001 shocked infall model and the most massive halos 
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Figure 4. The 2-8 keV (rest-frame) X-ray luminosity to SFR ra- 
tio at z = 6 as a function of host halo mass with r?o = 4. As in 
Fig- El solid (black), dotted (red), long-dashed (blue), dot-dashed 
(green), and short-dashed (cyan) curves denote results for LSW 
models with m = 0.2 and m = 1 and shocked infall models with 
(3 = 0.01, 0.1, and 0.001, respectively. The shaded, magenta re- 
gion marks the 2<r range of reasonable values for cx (representing 
emission from HMXBs) consistent with low-redshift observations 
of galaxies without AGN. In general, most of our angular momen- 
tum transport models produce a significant AGN X-ray excess 
intrinsically over the expected emission from HMXBs. 



in the m = 0.2 LSW model. These comparisons, of course, 
assume that the IMF and the resulting cx from HMXBs do 
not vary significantly at high redshift from their low-redshift 
values. If z > 6 galaxies have lower metallicities, the con- 
tribution to the X-ray emissi on from star formation could 
be higher than expected (e.g.. iBookbinder et al.|[l980l ; iDravl 
l200fj ; iLinden et al.l |2010| ; iKaaret et al.ll201lh . Additionally, 
the intrinsic X-ray emission from HMXBs in the star forma- 
tion disk and from a central AGN may be obscured in very 
different ways and to different degrees. 



6.2 X-rays from High-z LBGs 

Analyzing r e cent Chandra observations (|Xue et alj l201ll ). 
ICowie et alj (|2012l ) placed upper limits on the X-ray sig- 
na l from stacked popu l ations of LBGs at z « 6, 7, and 
8 (|Bouwens et al.ll2007l , fcoilbl ) in the soft, 0.5-2 keV (ap- 
proximately 3.5-14 ke V rest-frame) band. T hey also argued 
against the claims by iTreister et al.l (|201ll ) of a strong de- 
tection at z ~ 6. 

Using the probability distribution from equation l20l we 
randomly assign a halo mass to each galaxy in the LBG sam- 
ple. While the number of bright galaxies is not large enough 
to fully sample the halo mass distribution, their paucity also 
means that they contribute little to the average X-ray lumi- 
nosity despite each one's being individually much brighter 
than any of the more common UV-faint LBGs. We use these 



masses to calculate an X-ray luminosity (generated by our 
model AGN) for each LBG in the sample. We show the av- 
erage soft X-ray lu minosities in Figure [S] compared with the 
ICowie et all (2012) upper limits at z w 6-8. Also plotted 
for comparison is an estimate of the HMXB contribution to 
the X-ray luminosity from galaxies hosted in 10 10 Mq halos, 
where we have assumed an average SFR in these halos of 
0.2M Q /yr [(1 + z)/7] 2 - 5 , normalized to the UV LF at z = 6 
and with the same mass scaling as the cold-flow accretion 
rate. 

We find that the average X-r ay luminosity of z > 6 
LBGs is roughly consistent with the lCowie et al.1 l|2012h lim- 
its for all of our LSW models (left panel of Fig.0). Addition- 
ally, all nonlinear-shock models also achieve consistency if X- 
ray production is distributed through the plane of the AGN 
accretion disk so that the disk itself obscures the emission. 
By contrast, an unobscured, nonlinear accretion scenario 
only agrees with observations if j3 <C 0.01. The very rapid an- 
gular momentum transport rates of shock-induced accretion 
models produce, on average, about an order-of-magnitude 
more intrinsic X-ray emission than observed if /3 = 0.01 and 
about two orders-of-magnitude more if /3 = 0.1. We have 
only plotted results for r/o — 4 since setting 770 = 1 produces 
smaller changes than varying m, /3, or the amount of obscu- 
ration. Indeed, in the obscured models, the precise values of 
rjo and either m or /3, ultimately, have little effect on the 
resulting AGN luminosity. This is because changing param- 
eters to produce a higher BH accretion rate also tends to 
increase the absorbing column density. 



7 CONCLUSIONS 

We have extended the pressure-balanced ISM model of 
TQM05 to z > 6 using as few low-redshift or empirical pre- 
scriptions as possible. In this model, gas accretes onto the 
outer edge of the galactic disk and is transported towards 
the center. Along the way, star formation reduces the avail- 
able gas in an amount necessary to maintain hydrostatic 
equilibrium and marginal Toomre-stability. Deviating from 
previous studies, we assumed a dust-free ISM, included the 
disk self-gravity and the expulsion of a significant fraction 
of the accreted gas by winds. The galactic disk transitions 
smoothly into an accretion disk around a central BH provid- 
ing the energy to power an AGN. We calculate the physical 
and radiative properties of the disk as functions of radius 
by solving equation [16] under the constraints of equations [4] 
and [7] and assuming Q = 1 in the outer port ions of the disk. 

W e the n used t he UV LF analysis of iMufioz fc Loebl 
|201ll ) and iMunoj (|2012l ) to calculate the the distribu- 
tion of halo masses appropriate for z > 6 LBGs and cal- 
ibrate the wind mass-loading factor. A cold-flow accretion 
rate can readily reproduce the observati ons with r/o ~ 4 , 
consistent with numerical simulations l|Dave et al. I 120061 ; 
lOppenheimer fc Pavel I2OO8T ). However, where appropriate, 
we also explored the effect of assuming rjo = 1 where con- 
sistency with the UV LF is restored by a somewhat larger 
amount of dust extinction than currently indicated b y UV 
continuum slope measurements dBouwens et al.ll2007T l. For 
10 10 Mq halos at z = 6, these two values of 770 correspond to 
mass-loading factors of ?) w i n d ~ 8 and 2, respectively. While 
we do not investigate in detail a scenario in which major 
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Figure 5. Ave r age X- ray luminosity (0.5-2 keV observed frame; approximately 3.5-14 keV rest-frame) for the lBouwens et al.l i2007f ) and 
iBouwens et al.l l|2011bh LBG samples at z = 5.9, 6.8, and 8.0. Thick lines connecting filled symbols denote unobscured values and are 
contrasted with maximally-obscured luminosities shown in thin lines connecting open symbols. In the left panel, dotted (red, squares) 
and solid (black, circles) lines correspond to LSW models with m = 1 and m = 0.2, respectively. In the right panel, dot-dashed (green, 
squares), long-dashed (blue, circles), and short-dashed (cyan, triangles) denote nonlinear shocked infall models with /3 = 0.1, 0.1, and 
0.001, respectively. All of the calculations shown assume i)o = 4 as since effect setting r)o = 1 is much less than the difference between 
angular momentum transport mechanisms or between the maximally-obscured and unobscured cases. The shaded, magenta region marks 
the X-ray luminosities (corresponding to a 2<r range the distribution of ex) from a galaxy in a 10 10 M halo due to the HMXBs resulting 
from star formation. The horizontal, solid line shows the approximate ICowie et alj i2012h upper limits. 



mergers are responsible for z > 6 LBGs, we expect that 
such a mechanism would generate an unaccountably high 
star formation or BH accretion rate when compared with 
observations. 

A competition between the speed at which gas is trans- 
ported toward the center of the disk and the rate at which 
it is transformed into stars or ejected by winds determines 
how much gas remains to be accreted onto the central BH. 
Therefore, the distribution of star formation throughout the 
disk, the BH growth rate, and the resulting AGN luminosity 
critically depend on the angular momentum transport mech- 
anism in the disk. We compared the effect of three simple 
models: a local viscosity a-disk, infall mediated by a linear 
spiral wave (LSW), and nonlinear, shocked infall from orbit- 
crossings. We found that, while a-disks cannot account for 
a growth rate high enough to form a black hole commensu- 
rate with our assumed M-a relation, either an LSW model 
with m = 0.2 and r/o — 1 or a shocked infall model with 
P ~ 0.001 can build the correct size black hole by z ~ 6 with 
a duty cycle of unity (see Fig. [3]). We suggest that local, 
a-disk models and globally-operating, gravitational torque 
models (such as LSW and shocked infall) may represent dis- 
tinct growth modes for BHs in high-redshift LBGs. Further 
investigation is required to determine whether the rapid an- 
gular momentum transport mediated by dense clumps and 
disk instabilities (|Bournaud et al.l 120111 ) can, in a full cos- 
mological context, either stabilize or oscillate into a slower 
accretion mode dominated by local viscosity. 

Given the dearth of z > 6 AGN, we calculated the X- 



ray bolometric correction by extrapolating from low-redshift 
observations. However, the development of future samples 
will allow us to improve this procedure. While the X-ray 
emission generated in a-disk models is negligible due to 
minuscule BH accretion rates, AGN in both the LSW and 
shocked infall models generate more X-ray luminosity than 
do HMXMs. This is particularly true of shocked infall mod- 
els for all reasonable input parameters (the exception being 
for ft — 0.001), while an unusually high value of cx could 
allow HMXBs to compete with AGN in LSW models. How- 
ever, we also note that differences in the ways the AGN and 
the star-forming disk are obscured may complicate this com- 
parison. Future work, moreover, will show how this radiation 
feeds back into the IGM to modulate reionization. 

In considering the observable X-rays generated by our 
models, we considered two cases: unobscured emission-in 
which the emission is all produced at a radius of order the 
inner edge of the accretion disk and the disk is viewed face- 
on-and highly obscured emission-where X-ray production 
is distributed throughout the mid-plane of the accretion 
disk and radiation is spectrally-filtered through a disk scale 
height. Interestingly, we find that our more extreme mod- 
els are already ruled out by the current data; shocked in- 
fall models with f3 0.01 produce average intrinsic X-rays 
luminosities 1-2 orders-of-magnitude higher than observed. 
This leads to the conclusion that either angular momen- 
tum must be transported more slowly through the disk or 
the emission must be substantially obscured. While the re- 
maining models are all roughly consistent with the data, 
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none of the shocked infall models (obscured or unobscured) 
is more than an order-of-magnitude below the limits. The 
requirements of super-Eddington accretion and substantial 
obscuration may be circumstantial evidence arguing against 
non-linear, shocked infall in high-redshift LBGs. However, 
only moderately more sensitive observations are necessary 
to probe the more reasonable parameter- space. While each 
individual X-ray observation need not get deeper, stacking 
a larger population of high-redshift LBGs might bear fruit. 
S earch es to expand th e sample size of these objects are on- 
going (lYan et alj|20ld; iTrenti et al.ll201ll ; Ikoekemoer et all 
|20iJ :[g rogin et al.l l201ll ~ Focussing on lensed objects may 
also prove a way forward. However, it is clear that exploring 
high-redshift galaxies outside of the traditional UV bands 
will provide interesting opportunities to verify and inform 
models of the internal physics and formation of these sys- 
tems. 
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